Superconductivity Bose-Einstein condensation 



? 



S.K. Adhikari'', M. Casas'', A. Puente,^ A. Rigo,^ M. Fortes 
M.A. Soil's," M. de Llano,^ A.A. Valladares,'^ and O. Rojo'^ 



'^Institute de Fisica Teorica, Universidade Estadual Paulista, 01405-900 Sao Paulo, SP, Brazil 
^Deparment de Fisica, Universitat de les Illes Balears, 07071 Palma de Mallorca, Spain 
'^Instituto de Fisica, Universidad Nacional Autonoma de Mexico, 
Apdo. Postal 20-364, 01000 Mexico, DF, Mexico 

'^Instituto de Investigaciones en Materiales, Universidad Nacional Autonoma de Mexico, 
Apdo. Postal 70-360, 04510 Mexico, DF, Mexico 

^PESTIC, Secretan'a Academica & CINVESTAV, IPN, 04430 Mexico DF, Mexico 

Abstract 

Bose-Einstein condensation (BEC) in two dimensions (2D) (e.g., to describe the quasi-2D 
cuprates) is suggested as the possible mechanism widely believed to underlie superconduc- 
tivity in general. A crucial role is played by nonzero center-of-mass momentum Cooper pairs 
(CPs) usually neglected in BCS theory. Also vital is the unique linear dispersion relation ap- 
propriate to weakly-coupled "bosonic" CPs moving in the Fermi sea — rather than in vacuum 
where the dispersion would be quadratic but only for very strong coupling, and for which 
BEC is known to be impossible in 2D. 
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Bose-Einstein condensation (BEC) of Cooper pairs (CPs) leads to a phase transition even 
in 2D in any many-fermion system dynamically capable of forming CPs. This transition could 
be germane to superconductivity in the quasi-2D cuprates. In the weak coupling limit one 
finds a nearly linear dispersion relation for the CP that suggests very high, even diverging, 
critical temperatures Tc- On the other hand, in the strong coupling limit a nearly quadratic 
dispersion relation gives vanishingly small Tc's. For intermediate coupling one gets the finite 
Tc's appropriate for real quasi-2D superconductors. 

The single-Cooper pair problem may appear academic at first but has significant conse- 
quences. The familiar BEC formula for the transition temperature ~ S.Slh'^n'^J^ /mskB, 
with ub the number density of bosons of mass tjib and kB the Boltzmann constant, is a special 
case of the more general expression Q valid for any space dimensionality d > and any boson 
dispersion relation ek = Cg with s > and Cg a constant, given by 
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If fiB{T) is the boson chemical potential and e'^^^^^/'^^^ = z the fugacity, ga{z) = /I" ■ 

For z = 1 and a > 1 this is just C('^)) the Riemann Zeta- function of order a which is finite 
for cr > 1 and infinite for a = 1, while the series g'(r(l) diverges for all o" < 1. For s = 2, 
C2 = h?/2'mB, and since ^(3/2) ~ 2.612, this leads to the BEC Tc-formula cited above. Since 
5^/2(1) diverges for all d/2 < 1, Tc = for all d < 2. Eq. follows from the boson number 
equation 

N = No{T) + J2 [e{^^-^s(^)>/*=s^ - 1] (2) 

where No{T) is the number of bosons in the K = state. At T = Tc both No(Tc) and ^b{Tc) 
virtually vanish so that using ^jr. — > (L/27r)'^/ d'^k in (^) then yields (Q). 

The fact that a CP can have a linear, as opposed to the usual quadratic, dispersion relation 
was mentioned as far back as 1964 by Schrieffer p. 33. We have found it to be so for 
weak coupling; it becomes quadratic only for extremely strong coupling. In weak coupling 
ek = ^0 — — a{d)hvFK, where Ak is the (positive) binding energy of a CP of center-of- 
mass momenta (CMM) hK, vp = Tikp/rn is the Fermi velocity, ra the fermion effective mass, 
while a{d) = 2/7r and 1/2 in 2D and 3D, respectively |^. For linear dispersion s = 1, Ci = 
a{d)hvF, another special case of (|l]) is 
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Now Tc = for all d < 1 only — and > for all d > 1, which is precisely the range 
of dimensionalities for all known superconductors if one includes the quasi-lD organo-metallic 
Bechgaard salts Q. Using the interpolation a{d) = (7/2 - G/tt) + (S/vr - 13/4)d+ (3/4 - 2/7r)d^ 
which reduces to 1, 2/7r and 1/2 in ID, 2D and 3D, respectively. Fig. 1 graphs (|l|) and (^) vs 
d if one imagines all fermions in the many- fermion system paired into CPs of mass niB = 2m. 
The particle number density of the original fermions is n = 4/2'^-27r^/2^r(tf/2) = 2nB. Here 
kBTp = Ef = h'^kl/2m. These curves are an upper bound to a more realistic Tc where only a 
fraction of all fermions are actually bound. 

The above holds for "unbreakable" CPs, i.e., bosons existing for all CMM so that the integrals 
over momenta go from to 00. However, CPs actually break up |^ into two "unpaired fermions" 
above a certain CMM Kq. This increases in value as coupling is increased and tends to infinity 
(as will be seen below) for infinite coupling. The CP breakup Kq, defined by Akq = 0, for very 
weak coupling is just Kq ~ AQ/a{d)hvF ^ as Aq — 0. Eq. (|2|) with Kq rather than cxd as an 
upper limit eventually yields 

^ ^ {d - l)[27ra{d)hvF]''nB _^ ^ 
A^ksA^^-' Ao-o ' 

where = 27r or 47r/3 in 2D or 3D, respectively. Hence, since in that limit any remaining 
CPs are condensed in the K = state at all temperatures, the critical temperature Tc must be 
infinite. Note that Tc = explicitly for d = 1. Indeed, if the unpaired fermions resulting from 
broken CPs are taken into account in a more realistic binary boson-fermion mixture model this 
artificially infinite Tc is "tamed" down Q to a finite Tc — which remarkably enough practically 
coincides with that of the pure unbreakable boson gas (H). Detailed calculations using the 
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Figure 1: Upper bounds to BEC critical temperatures Tc as function of system 
dimensionality d, as explained in text. Shaded areas refer to empirical data Q while 
BCS "phonon barrier" refers to BCS ~ 1.138^ exp(-l/A) with Gd the Debye 
temperature ~ 300 K and A = N{0)V = 1/2 where A^(0) is the density of states at 
the Fermi surface and V the attractive BCS coupling constant. 



BCS model interaction gives a Tc still almost an order of magnitude too large compared with 
the empirical range of Tc ^ (0.01 — 0.1)rp for cuprates — even for a moderate coupling of 
A = 1/2. However, the calculated Tc of about 800 K is about 17 times larger than that of the 
BCS theory Tc ~ 1.13©d exp(— 1/A) of about 46 K. For cuprates d ~ 2.03 has been suggested 

as more realistic as it accounts for inter-layer couplings, but results would be very similar to 
those discussed here for d = 2. 

A far more general interfermion interaction is the S'-wave attractive separable potential 0] 
whose double Fourier transform is 

Vpq = - {vo/L'^)9p9q- (5) 

Here L is the size of the "box" confining the many-fermion system, > is the interaction 
strength and = {l+p"^ /Pq)~^^^ where po is the inverse range of the potential. Hence, e.g., po — >• 
oo implies = 1 and corresponds to the contact or delta potential V{r) = —voS{r), andpo = kp 
to a range of order of the average interfermion spacing, etc. Such an interaction model may 
mimic a wide variety of possible dynamical mechanisms in superconductors: a force mediated 
by phonons, or plasmons, or excitons, or magnons, etc., or even a purely electronic interaction. 
In the first instance mentioned one may have, e.g., a (possibly singular) coulombic interfermion 
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repulsion surrounded by a longer-ranged electron-phonon attraction. The former is central to 
the problem at hand; the latter (whether or not phonon-based) is indispensable to create the 
CPs detected experimentally |8| in both conventional as well as cuprate superconductors. In 
vacuum, a two-body bound state with (positive) binding energy i?2 ^ is possible for any 2D 
potential with sufficient attraction. The vacuum t-matrix Q associated with the Lippmann- 
Schwinger equation then develops a pole at an energy E = —B2 so that for any dimension d @ 
leads to 

Td 2 

_ = V ^ f6) 

In ID the integral on the rhs is exact and gives the well-known binding energy B2 = mv^ / 4^^ of 
the single bound state, but in 2D that integral diverges logarithmically for large k. On the other 
hand, the CP equation for two fermions above the Fermi surface (of energy Ep) with momenta 
wavevectors ki and k2 (and CMM wavevector K = ki + k2) is given by 

V' ^-^ = — (7) 

Y n'^k^/m + /:\K-'2.EF + n'^K'^/^m Vo 

where k = ^(ki — k2) is the relative momentum wavevector, 2Ep — Ak the total pair energy, 
Ak > as before, and the prime on the summation implies restriction to states above the Fermi 
surface, viz., |kzbK/2| > kp. In principle, (0) can be solved numerically for Ak as a function of 
coupling strength t^o for all < f;o < +00. More conveniently, however, for an interaction of the 
form dsj) one can solve for Aq (analytically if gk = 1) and for Ak numerically, both as functions 
of B2 alone, after eliminating the interaction strength vq by combining and (^). This yields a 
renormalized CP equation em.\Aay\ng < B2/EP < 00 [^] as dimensionless interaction coupling 
parameter instead of the (possibly singular) potential strength term I/vq is now eliminated. A 
much greater generality for the equation with gj. = 1, viz., for any interfermion interaction with 
binding energy B2, is conceivable. In 2D one has the surprising result (traceable to the fact that 
unlike ID or 3D the fermionic density of states is a constant independent of kp) that Aq = B2 
for all coupling, at least for an attractive delta interaction. Fig. 2 shows exact numerical results 
for the nonnegative CP excitation energy sk = Aq — A^ mentioned earlier, for zero range (full 
curves). For very weak coupling the exact curves are virtually linear, i.e., £k 2hvpK/ir, and 
for very large coupling tend aymptotically to the exact quadratic fi'^K'^ /2{2m) (short-dashed 
curves). The long-dashed loops (reminiscent of the "roton" part of the excitation spectrum of 
liquid ^He) emerge for finite range at sufficiently strong coupling. 

When the zero-range results like those in Fig. 2 are used in (|2|) assuming ns = 'n/2 and 
niB = 2m, one gets Fig. 3 for the BEC Tc's of a pure gas of breakable or unbreakable CPs. 
Consequently, for a substantial range of intermediate couplings in 2D finite BEC Tc's are possible 
that lie within the empirical range Q of cuprate Tc's even for a simple pure boson gas of CPs. 
Significantly, Tc is no longer zero as would be predicted in a BEC picture by a quadratic relation 
appropriate for "local-pair" CPs in vacuum — a result that convinced many that BEC is irrelevant 
for quasi-2D cuprate superconductors. More accurate BEC Tc's should include refinements such 
as allowing for unpaired fermions in a more realistic binary boson-fermion mixture model, 
CP-fermion interactions, non-S-wave interactions, etc. 

Finally, CPs are here considered "bosonic" even though they do not obey (Ref. p. 38) Bose 
commutation relations. This is because for a given K they have indefinite occupation number 
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Figure 2: Dimensionless CP excitation energy £k/^o vs dimensionless CMM 
wavenumber K/kp, as explained in text. Full curves are exact zero-range results; 
short-dashed the quadratic approximation; and long-dashed the exact finite-range 
result with po = kp- Dots mark CP breakup points. 



as (in the thermodynamic limit) there are an indefinitely large number of allowed (relative 
wavevector) k values corresponding to an indefinitely large number of possible pairs of vectors 
ki and k.2. Hence, for any coupling and thus any degree of overlap between them, CPs do in 
fact obey the BE distribution ^ from which BEC is determined. 

To conclude, the single CP problem with non-zero CMM (usually neglected in BCS theory) 
is illustrated in 2D for all coupling (and/or normal-state fermion density). BEC then suggests 
itself as a possible mechanism for superconductivity — in fact, for all observed superconductors 
from the quasi-lD organo-metallics to the more familiar 3D ones. But one must employ the 
correct dispersion relation in the CMM of the CPs. This relation is complicated but obtainable 
numerically for arbitrary coupling, and is purely linear for very weak and purely quadratic for 
very strong coupling. 

Partial support from UNAM-DGAPA-PAPIIT (Mexico) # IN102198, CONACyT (Mexico) 
# 27828 E, DOES (Spain) # PB95-0492 and FAPESP (Brazil) is acknowledged. 
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Figure 3: Critical BEC Tc (in units of Tp) calculated from (^) for a pure gas of break- 
able (i.e., with finite limit Kq) CPs using the exact numerical dispersion illustrated 
in Fig. 2, vs all coupling < B2/EF < 00. The unbreakable (i.e., Kq = 00) 
case corresponds to the interpolation ansatz ek = -fiypKll — tanh(i?2/-E'F)] + 
^Unh{B2/EF). 
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